Confinement effects on the stimulated dissociation of molecular BECs 
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We show that a molecular BEC in a trap is stabilized against stimulated dissociation if the trap 
size is smaller than the resonance healing length (h 2 /2mg^/n) 1 ^ 2 . The condensate shape determines 
the critical atom-molecule coupling frequency. We discuss an experiment for triggering dissociation 
by a sudden change of coupling or trap parameters. This effect demonstrates one of the unique 
collective features of 'superchemistry' in that the yield of a chemical reaction depends critically on 
the size and shape of the reaction vessel. 



Atom-molecule conversion in quantum bose gases is 
equivalent to three-wave mixing in nonlinear optics 
[2 S B EL H, @1- Association is analogous to second- 
harmonic generation and dissociatiocorresponds to para- 
metric downconversion. The latter process exhibits a 
modulational instability, resulting in from exponential 
amplification of spontaneously emitted photon-pairs 0] • 
The existence of similar parametric gain in molecular 
dissociation was noted early on in the study of atom- 
molecule condensates 0, 0, 0, 0]. This dynamical in- 
stability generates significant deviations from the Gross- 
Pitaevskii (GP) classical field theory [1, [n|. The rapid 
growth of correlations can be used for the generation of 



pair-correlated or number-squeezed atomic beams jllj in 
analogy to the use of parametric downconversion to gen- 
erate squeezed light [12| . 

Exponential gain in molecular BEC dissociation is 
readily ca ptu red within an undepleted pump approxima- 
tion (3. lil. Illj. The atom-molecule Hamiltonian is, 
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where if) and (j) are atomic and molecular bose field op- 
erators, m is the atomic mass, A is the detuning from 
atom-molecule resonance and V(r) = muiQ 2 r 2 /2 is an 
isotropic magnetic or optical trap potential with fre- 
quency ljq. Atoms and molecules see essentially the same 
trap frequency since molecules have twice the mass and 
twice the magnetic moment or optical polarizability [l3| . 
Atom-molecule coup ling of strength g, is achieved by a 
Feshbach resonance 14[ or by an optical stimulated Ra- 
man transition [IB}. The molecule- molecule interaction 
strength is U m and interactions between atoms are ne- 
glected since we are interested in the initial dissociation 
dynamics when atomic densities are small. The Heisen- 
berg equations of motion for the field operators read, 
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The undepleted pump approximation is a linearization of 
Eq. ([2]) about a stationary, classical molecular field <fi(r) 
obeying the GP equation 

H 2 

-— V 2 (j ) +[2V{Y) + U m \^\ 2 ]<t> = ^ . (4) 
4m 

Replacing <j){v,t) — > (f>(r)e~' lfJ/t / h , and rotating ^ = 
^ e v*/2^ = ^t e -v«/2fi ( so that the atomic density 

remains n a = ffi'ip) = we obtain the coupled 

linear equations for the atomic field operators *P and 
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where S = A — yu/2 is an effective detuning including 
the molecular chemical potential consisting of zero-point 
energy and of a molecular interaction shift. The molecu- 
lar condensate stability is determined by the eigenvalues 
A of the set ©-([B]). For a uniform gas with V(r) = 0, 
<p = ypri, and ji = U m n (n being the molecular den- 
sity), the system is diagonalized in momentum space 
i/j(r, t) = J^q Oq(i) exp(iq • r), to give the dispersion, 



h\ q = ^(e q + 5) 2 -ng 2 
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where e q = (hq) 2 /2m is the free-particle dispersion. Thus 
for a uniform gas with 5 < \/n\g\, there are unstable 
resonant modes for which e q + S < g\fn, resulting in 
complex characteristic frequencies and exponential gain. 
The dynamics of amplified modes is described by well- 
known solutions, given in terms of the hyperbolic func- 
tions cosh(Agi) and sinh(A 9 t), which can be found in Rcfs. 
[1, EI , El I ■ For 5 = 0, the unstable modes are simply low- 
energy, long-wavelength excitations. 

To gain insight on nonuniform molecular BEC stabil- 
ity, we c omp are Eqs.©-© to the standard Bogoliubov 
analysis [l6| of an attractively interacting BEC, with the 
uniform-gas dispersion X q = \J £q{t q + 2nU), where U is 
the interaction strength. The difference between the Bo- 
goliubov dispersion relation and Eq. ([7|) results from the 
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FIG. 1: (color online) Real and imaginary parts of character- 
istic dissociation frequencies as a function of atom-molecule 
coupling strength, for a harmonic trap with frequency u>q. 
The coupling is assumed to be resonant, i.e. 5 = and the 
molecular condensate is noninteracting, i.e. U m = . Com- 
plex frequencies only appear when the coupling strength g\fn 
is larger than the trap level spacing huio 



existence of additional mean-field diagonal terms in the 
Bogoliubov equations, corresponding to different pairings 
of the Uij^^ij) interaction term. For attractive interac- 
tion U < 0, the long-wavelength modes with e q < \2nU\ 
are unstable, leading to BEC collapse [l7j]. For a BEC 
confined to a trap however, collapse is prevented pro- 
vided that the size of the trap is short with respect to 
the healing length £ = y / h 2 /(2mnll). Since the unsta- 
ble modes are those excitations with wavelength longer 
than £, they are precluded from tighter traps. Therefore, 
stable condensates of attractively interacting atoms exist 
up to a critical density for a given trap-frequency [lij ]. 

A similar instability threshold exists for molecular 
BEC dissociation. Even for (5 = 0, there will be no 
unstable modes towards dissociation provided that the 
trap size lo = ^/h/JmLUo) is smaller than the 'resonance 
healing length' ( = y h 2 / (mgn 1 / 2 ), which is the charac- 
teristic lengthscale of Eqs. ([5|) and ([6]), resulting from 
balancing kinetic and resonance energies. Thus for suf- 
ficiently tight confinement, atomic zero-point motion in- 
hibits dissociation. Remarkably, the outcome of a coher- 
ent chemical process depends due to its collective nature, 
not only on the statistics of the constituent atoms, but 
also on the size and geometry of the 'container' in which 
it is carried out. 

In order to ascertain that confinement can stabilize a 
molecular BEC against stimulated dissociation, we have 
calculated the characteristic frequencies of Eqs. ([5])- 
(J6j) with a finite-size undepleted pump <f> obtained from 
Eq. Q, by expanding 0(r) = EjVjXjWj ^( r > 4 ) = 
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FIG. 2: (color online) Total atomic fraction obtained from 
a mean-field calculation, triggered by Gaussian noise, as a 
function of rescaled time r = ujot for g^/n = 0.9huJo (dashed 
red line) and for g^/n — fiuot (solid blue line), with 8 — U m = 
0. The corresponding atomic density profiles at the final time 
r = 30 are shown in the inset. 



Ej % (*)Xj ( r ) m the basis of harmonic-oscillator solu- 
tions Xj( r )i an d diagonalizing the resulting linear set of 
equations for the dynamics of dj(t) and a At), 
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with a = (a ,ai, ...,a,j, ...), &< = (%,&[, ...,aL ...), fiy = 
(j + 1/2), $ y = (ff/RJEjWteilXi)- The re- 
suits were compared to a discrete Fourier transform di- 
agonalization, giving precisely the same eigenfrequen- 
cies. Eigenvalues of a one-dimensional calculation with 
S = U m = 0, are plotted in Fig. [T] as a function of 
the atom-molecule conversion frequency g^/n/h. Com- 
plex frequencies, corresponding to an initially exponen- 
tial molecular gain, are only obtained if the coupling 
frequency is larger than the trap frequency. Thus, for 
g^/n/fiWQ < a where a is a geometrical factor of order 
one, the molecular condensate is stabilized by atomic zero 
point motion. 

In Fig. [2] we plot the results of an approximate 
solution of the dynamical equations (JH)-©, including 
pump depiction. The equations were solved in a mean- 
field approach, simulating initial spontaneous dissocia- 
tion by a complex Gaussian white noise with zero mean, 
corresponding to atomic quantum fluctuations in the 
Wigner representation [rJ [2(j. While this approach 
does not fully capture the quantum dynamics in the way 
full stochastic calculations do [ll|, it suffices to depict 
the dominant quantum effect of amplification of spon- 
taneously emitted atom pairs in the vicinity of the in- 
stability 0, d]. Since we are interested here in the 
mere presence of an instability, our calculation should 
indicate whether molecular gain exists even if it does 
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FIG. 3: (color online) Critical coupling strength ^/ng c as a 
function of the molecule-molecule interaction frequency nU m . 
Magnetic or optical atom-molecule detunings A were assumed 
to compensate the molecular chemical potential, so that 5 = 
for each interaction frequency. In this way, true shape effects 
are distinguished from mean-field shifts originating from the 
variation in the molecular chemical potential \i. 



fail to give the precise onset time. Evidently, Fig. [2] 
demonstrates the existence of an amplification threshold 
g c ps 0.91frw / ' \fn. For g < g c all frequencies are real 
and there is no observed dynamical gain, whereas for 
g just above the critical value, molecular amplification 
takes place. 

The precise critical ratio between coupling strength 
and trap frequency depends upon the shape of the molec- 
ular condensate. In Fig [3] we plot the numerically- 
obtained value of a as a function of the molecule- molecule 
interaction strength nU m , in a one-dimensional calcu- 
lation. The condensate profile becomes broader with 
increasing repulsive interaction (see inset), leading to 
higher values of g c for a given uj and n. Since the cal- 
culation is in ID, stable molecular condensates with self- 
trapped solitary profiles, significantly narrower than the 
trap size, do exist also for attractive interaction. The 
lowest value of a is obtained at the transition to self- 
trapping at nU m ~ —hujQ. 

The coupling energy scales as g^/n = g(N / Kdlo 1 ) 1 ^ 2 
where d is the number of dimensions and Ad is a d- 
dependent geometrical factor. Since the spacing between 
trap levels is Hujo — h 2 /mlo 2 , we obtain that at the crit- 
ical point ^N/k d lv 2 ~ d/2 ps h 2 /{mg). Thus for a three- 
dimensional system d = 3, A3 ps 8n^/2/3 and a fixed 
trap size, the critical number of particles in the molec- 
ular condensates is Nc ~ A3 C^^j /lo- This critical 
number should be compared with the critical number of 
particles for the collapse of a 3D BEC with attractive 
interaction of strength U = Anh 2 a s /m < (a s being the 
atomic s-wave scattering length), obtained from balanc- 
ing 2Un = 2UN/(A 3 l 3 ) with huj = h 2 /ml 2 , to give 



the critical number of particles for collapse, larger molec- 
ular condensates will be stabilized against stimulated dis- 
sociation for tighter traps, because the increase in zero 
point motion overcomes the effect of increasing density. 

Most atom-molecule systems coupled by a magneti- 
cally tunable Fcshbach resonance, are well within the 
strong-coupling, nonlinear domain, where the character- 
istic atom-molecule conversion frequency is much larger 
than the trap frequency. For example for the 23 Na Fesh- 
bach resonance at 907 G 21 ,23], the coupling strength 
is given by g = \/nU a |23l. 24] with U a = A-Kh 2 a s /m, 
k/(2ttK) ps 4.6 MHz, and a s ps 60 Bohr. These val- 
ues give g/2Trh ps 9 x 10~ 3 Hz cm 3 ' 2 . Given a density 
of n = 10 15 cm~ 3 , we have g^/n/2irh ss 284 kHz, far 
above contemporary trap frequencies (we note that com- 
parison of this Fcshbach frequency to the atomic interac- 
tion frequency U a n/2iiTi ps 17.5 kHz, justifies the neglect 
of atom-atom interactions in our calculations, even for 
substantial atomic populations and high densities. For 
n = 10 13 cm~ 3 , gy/n/2irh drops to 28.4 kHz, whereas 
U a n/2irh is only 175 Hz). For cj /2tt = 200Hz we ob- 
tain that N c Ri 2 x 10 -2 , so that stimulated dissociation 
will be observed for any number of molecules in the con- 
densate. It may be that a sufficiently weak Feshbach 
resonance can be found amongst the many resonances of 
87 Rb 0. 

An experimental demonstration of collective dissoci- 
ation thresholds, can however be facilitated in setups 
where atoms and molecules are opti cally coupled via a 
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stimulated Raman transition [26 
this case, g = £l p £l s /2A2 P is an effective two-photon Rabi 
frequency. The one-photon Rabi frequencies f2 SiP for the 
pump (inducing a bound-bound transition to an interme- 
diate molecular state) and the Stokes (dissociating the 
intermediate state via a bound- free transition) lasers, re- 
spectively, are products of laser intensities by overlap in- 
tegrals consisting of electronic transition dipole moments 
and Franck-Condon factors, and A2 P is the two-photon 
detuning from the intermediate bound state. Thus, the 
effective coupling strength can be controlled with great 
precision [28] through the adjustment of laser parame- 
ters. Most app ealing are systems of molecules in opti- 
cal lattices 27, 28, 3fl 3l|], providing optical coupling to 
deeply bound internal molecular states, as well as tight 
trap frequencies of luq/2-k — 10 — 100 kHz. Currently, 
these experiments seek to avoid collective effects in as- 
sociation by operating in the Mott-insulator regime with 
unit occupation numbers. However, provided that molec- 
ular condensates containing roug hly 10 2 - 10 3 particles 
per site be formed, their dissociation can demonstrate 
the expected amplification threshold. 

Given 10 3 87 Rb2 molecules in a 1 kHz trap, the critical 
value g c = ti 2 / (my NI0/A3) of the interaction strength 
turns out to be approximately 2t:Hx 2.2 x 10 -5 Hz cm 3 / 2 . 
The corresponding average atomic density is of order 
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10 14 cm" 3 and a s ps 100 bohr for 87 Rb. Therefore, 
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FIG. 4: (color online) Dynamics of dissociation following a 
sudden switch in g from 0.9?kJo/ \/n to huio/y/n (a-b) or in u>o 
from l.lg^/n/h to g^/n/h (c-d). Dashed lines mark the time 
of the switch at u>ot = 4n. Atomic density distributions are 
shown in (a) and (c) using a logarithmic colormap, whereas 
the molecular condensate density is plotted in (b) and (d) on a 
linear colormap. Molecular interaction and effective deuning 
are U m =5 = 0. 

Un ft 75 — 750 Hz for a fully dissociated gas, and atomic 
interactions are initially negligible with respect to the 
trap and coupling frequencies. The proposed experiment 
will involve the preparation of a molecular BEC in a trap 
smaller than the resonance healing length, thus arresting 
dissociation even in the presence of coupling. A sudden 
change either in g or in ujq (which in a lattice can be 
achieved via switching the lattice wavelength or depth) 
will trigger the stimulated dissociation of the BEC. Both 
scenarios are depicted in Fig. [4] for a ID calculation. 
In Figs. [4jt and [4j), atomic and molecular densities are 
shown as a function of time for a <?-switch experiment. 
Evidently, dissociation is triggered by the change in the 
coupling frequency. Similarly, a sudden reduction in trap 
frequency (Figs. SJ:,d) can lead to stimulated dissociation 
of the breathing molecular condensate. When the simula- 
tion is carried over the same timescale without changing 
g or ljq, no dissociation is observed. 

In conclusion, progress in tight confinement of molec- 
ular quantum gases and control over atom-molecule cou- 
pling, brings trap frequencies to the vicinity of atom- 
molecule conversion frequencies. Molecular BECs can 
be stabilized against Bosc-stimulatcd dissociation, if con- 
fined to a trap smaller than the resonance healing length 
£. Critical values depend on condensate shape and trap 
geometry. Experimentally, dissociation can be triggered 
by a sudden change in trap parameters. Such an exper- 
iment requires that (a) £ ~ lo to arrest dissociation, (b) 
uIq 3 > 1 so that collective effects are significant, and 
preferably (c) £ 3> C so that inelastic collisions between 
atoms do not interfere with association. Future work will 
explore pattern formation in 2D and 3D tight-trapping 



experiments. 
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